Two spin measurements in exchange interaction quantum computers 
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We propose and analyse a method for single shot measurement of the total spin of a two electron 
system in a coupled quantum dot or donor impurity structure, which can be used for readout in a 
quantum computer. The spin can be inferred by observing spin dependent fluctuations of charge 
between the two sites, via a nearby electrometer. Realistic experimental parameters indicate that 
the fidelity of the measurement can be larger than 0.9999 with existing or near-future technology. 

PACS numbers: 03.67.Lx, 73. 21. La, 72.25.Rb, 03.65.Ta 
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Semiconductor technology is rapidly reaching the level 
where quantum systems comprising of one or two elec- 
tron spins can be confined and coherently manipulated 
within a single nanostructure Q, • These experiments 
are important from a fundamental point of view, and will 
also pave the way to new applications in quantum infor- 
mation processing (QIP) 0. [J 13 A key requirement 
for these QIP schemes is the realization of a single shot 
readout technique, whereby the spin state of a one- or 
two-electron system can be determined in a single mea- 
surement run. Such measurements are important both 
for the ongoing experimental development of spin qubit 
systems (e.g. for characterizing qubit parameters, and 
studying the physics of decoherence processes) and, ulti- 
mately, for developing scalable QIP architectures. 

Recently, elegant experiments have demonstrated sin- 
gle shot readout of single electron spins,yia optical tech- 
niques in a defect center in diamond and using an 
all-electrical technique involving spin to charge conver- 
sion in a quantum dot, together with charge detection 
with a nearby quantum point contact (QPC) |2|. A num- 
ber of other methods based on spin-to-charge conversion 
have also been proposed 0, 0, 0| • Of particular inter- 
est are readout schemes in donor impurity implementa- 
tions of QIP 0, 0| , in which an electric field is applied 
to induce spin-dependent polarization of a two-electron 
double-donor system. It is believed that these propos- 
als may fail due to the short lifetime of the quasi-bound 
two electron system under the large E field required to 
observe a significant polarization j^, llflj . In this paper 
we discuss an alternative scheme, whereby the total spin 
of a two electron system can be inferred by observing 
spin dependent fluctuations of charge between two tun- 
nel coupled quantum dots (CQD) |22j, even in the ab- 
sence of an external electric field. Our scheme therefore 
avoids the lifetime issue in donor impurity systems, and 
furthermore may be easily integrated into proposed QIP 
architectures, as it does not require magnetic tunnel bar- 
riers, tunnel coupling to electron reservoirs, or external 
R.F. fields. The scheme can be used for read out of a sin- 
gle spin qubit (together with an ancilla spin), read out 
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FIG. 1: A schematic of the singlet-triplet readout scheme. 
Configurations (a) and (c) are energetically forbidden for 
triplet states, and thus monitoring fluctuations in the QPC 
current allows a measurement in the singlet-triplet basis. 



of an encoded qubit [llj . or, together with appropriate 
single spin operations, as a two-qubit readout device. 

Spin dependent charge fluctuations can be observed by 
continuously measuring the CQD system with an elec- 
trometer adjacent to one of the dots (see Fig. For 
the purposes of this work, we model the electrometer as a 
QPC, although it might equally be a single electron tran- 
sistor. Note that the QPC has a dual role in this scheme: 
it acts both as a noise source which can induce inelas- 
tic transitions in the CQD system, and as a detector to 
observe these transitions. 

In what follows, we first describe a model for the cou- 
pled CQD-QPC system, and subsequently derive master 
equations for both the unconditional and conditional dy- 
namics of the CQD system. We use these results to sim- 
ulate individual runs of the measurement scheme, and 
to characterize the detector output for the different mea- 
surement outcomes. We show that observing the detector 
output leads to a quantum measurement in the singlet- 
triplet basis: for a singlet state (total spin S = 0) one 
observes fluctuations in the output current, while for a 
triplet state (S = 1), these fluctuations are energetically 
suppressed, owing to the Pauli principle. Finally we de- 
termine the single-shot singlet-triplet measurement time 
and fidelity for realistic experimental parameters. 

Model. A number of authors have considered the prob- 
lem of a continuously observed single charge in a coupled 
dot system 0, [lj, [lj, [lj, [l5| . Here, we adopt the quan- 
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turn trajectories approach of 0, which takes into 
account the role of energy exchange mechanisms between 
the observed system and a QPC detector at finite voltage 
bias. We describe the two-electron coupled-dot system by 
a two site Hubbard Hamiltonian of the form 



dynamics for the CQD system is given by 



H = hA(t) 



where a^J is the fermionic annihilation (creation) oper- 
ator for an electron on site i with spin a, <tii a — a ia a,i ai 
A(t) is the tunnelling amplitude between sites, and U is 
the on-site Coulomb energy. We allow A(t) to be time 
dependent, since this can lead to shorter measurement 
times (see below), although in what follows, we frequently 
suppress this time dependence for clarity. H is spanned 
by the four singly occupied states \0\02) (where Oi = f, J, 
denotes the spin on dot i) and two doubly occupied states 
\di) = at^ajjO) (i = 1,2). The eigenstates of H are as 
follows. The singly occupied triplet states \t + ) = |||), 
|*d> = 2-5(|TI) + |lt)) and |t_) = HI) form a degenerate 
subspace with eigenvalue 0. The remaining states (with 
total spin S = 0) are |s ) = 2~* cos(0/2)(|U) - |||)) + 
2-ism(9/2)(\d 1 ) + \d 2 )), \ Sl ) = 2"* sin(0/2)(|U) - 
lit)) + 2-3 cos(0/2)(|di) + |d 2 )), and |s 2 ) = 2-^{\d x ) - 
I gJ 2 ) ) , where 8 = tan _1 (4A/J7). These states have eigen- 
values — J, U + J and U, respectively, where J(i) = 
(v/C/ 2 + 16A 2 (i) - U)/2 w AA 2 (t)/U is the exchange 
splitting between the lowest energy singlet and triplet 
states. J(t) can be varied over many orders of magni- 
tude by varying a surface gate voltage |3J, LLa] • 

The entire system (i.e. CQD coupled to a QPC) has a 
Hamiltonian H tot = H + H tU n + -fficads where 

fftun = h ( T kq + XkqTii) a\ k<y a Rqa + h.c. , (2) 

k,q,a 

H icads = h^2 u Lk a[ kiT a Lka +hJ2 UR q (J Rqa a Rqa . (3) 



Y (a\ a a 2(7 + al a a la ) + hU ^ * 

7=1,1 »=1,2 
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Here, a Lkr7 (a R ) are fermionic annihilation operators 
for electrons in the fcth (qth) mode on the left (right) 
side of the barrier with spin er and angular frequency 
^Lk i^B.q)- Tunnelling between the left and right baths is 
described by the factor T kq +\ kq ni, where n\ = ni-f +nxi 
is the total occupancy of dot '1'. T kq « T av and \ kq w 
Xav are assumed to vary slowly over the energy range 
where tunnelling is allowed. 

Unconditional master equation. We now derive an un- 
conditional master equation (UME) for the evolution of 
the reduced density matrix of the CQD system, p(t). 
The UME is obtained by first transforming to an inter- 
action picture in which the dynamics are governed by 
Hi(t) = W(t)H tun U(t), with U(t) = e -i(ff+iW s )t. l n 
the Born and Markov approximations |l7j . the resulting 
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dt' [Hi (t) , [Hi (t' ) , pi(t)®p L ®p R ]]}. 



(4) 

Here, pi(t) = (t) p(t)U (t) , and pL and p R are 
the lead density matrices, given by Fermi-Dirac 
distributions with chemical potentials pl and pu. 
The interaction picture tunnelling Hamiltonian may 
be written explicitly as Hj(t) = J2k. q .A T k q + 
Xkqni(t)]e t( > u>Lk < T ~ ujRq ''^ t a^ Lkij aR q(T + h.c, where the time 
dependence of n±(t) is given by n\{t) = ^ — 



,-i(J+U)t 



sin(0/2)|s o )(s 2 | 



cos(9/2)\ Sl )( S2 \ 



Note that Eq. J3J corresponds exactly to the first non- 
vanishing term in a perturbative expansion of pi(t) in 
terms of Hi(t) |l7j . 

We now make two further controlled approximations. 
First, we make a rotating wave approximation, setting 
very rapidly oscillating factors e lait — > 0, where on = 
{U,J + U}. This is valid for a, ^> v 2 V ', where v = 
y/^gLdRXnv is a dimensionless coupling constant, is 
the density of states in lead i, and V = (pl — PrJ/Jt- 
Second, we neglect small terms of order v 2 \J JV , which 
is valid when J -C V. At low temperatures (fc^T <C ?iV), 
and taking V > J + U, the UME is 



p{t) = -- [H,p(t)} +J2v[c t ]p(t) = Cp(t) , (5) 



where T>[c]p = J[c]p — A\c\p, J[c]p — cpc\ A[c]p 
(c'cp + pc'c)/2. The Lindblad operators are given by 



ci = vyJV - {J + Uj sin -\s 2 )(s \, 



c 2 = vs/V+(J + U) sin-|so)(s 2 |, 



C3 



T + f - z/cos -(|si)(s 2 | + h.c.) 



(6) 
(7) 
(8) 



where T = v^ffifffl^av The c^'s can be associated with 
different types of tunnelling process in the QPC. ci^ cor- 
respond to inelastic transitions, in which electrons tun- 
nelling through the QPC exchange energy with the CQD 
system, and are accompanied by transitions between the 
low-energy, singly occupied state, and the high energy, 
doubly occupied state. C3 corresponds to a quasi-elastic 
transition, in which electrons tunnel through the QPC 
without changing energy. 

Conditional dynamics. Individual measurement runs 
can be simulated using a conditional master equation 
(CME), which describes the evolution of the CQD sys- 
tem conditioned by the observed detector output, and 
also permits calculation of the current power spectrum. 
To derive the CME, we use an explicit model of the mea- 
surement process in terms of projective measurements of 
the number of electrons that have tunnelled through the 
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FIG. 2: A sample evolution of the conditional density matrix, 
p c , during a single measurement run. (a) Curves i = {(i) 
. . . (iv)} denote the probabilities P(i,t) = tr[|0i) {<f>i\p c (t)] for 
the system to be in the state \4>i){4>i\ at time t, where \4>i) = 
{|TD, UT), \di), \d%)}. (b) Corresponding probabilities for 
the system to be in a triplet state (curve (i)) or singlet state 
(curve (ii)). (c) Corresponding QPC current, as a histogram 
of tunnelling events. The parameters used in this simulation 
are A = 0.05(7 (we keep A(t) constant for this simulation), 
V = 5J7, T — \/03 « 0.54, v = -0.15T « -0.082. The initial 
state is taken to be p c = |TI)(T1|- In this measurement run, 
the final state of the system is a singlet state. 

PC, similar to that presented in [H [T^]. The CME is 
found to be 



dp c (t) = dN c (t) 



Plc(t) 



Pc{t) 



dt 



■ ^ ACn]Pc{t) + Plc(t) + l[H qbl Pc {t)] 



(9) 



Here, dN c (t) is a stochastic point process, correspond- 
ing to the number of electrons (0 or 1) which tun- 
nel through the QPC in an interval dt, P\ c {t)dt = 
t f {J2 n J[ c n\Pc{t)} dt is the corresponding probability of 
observing such a tunnelling event. Note that, in deriving 
Eqs. (J5J and (JSJ, we have used a course graining in time, 
which means that these equations are valid for timescales 
longer than U -1 . In practice U > 10 12 s _1 , and so this 
is not a significant restriction. 

Two sample solutions of Eq. @ are shown in Figs. [21 
and|3| In Fig. Ufa) the system initially undergoes almost 
coherent oscillations between the singly occupied states 
III) (Til an d lit) (IT I at the exchange frequency J. Af- 
ter a certain time the dynamics undergoes a qualitative 
change: the system suddenly jumps into the doubly oc- 
cupied state and thereafter undergoes stochastic 
jumps between states with distinct charge configurations. 
This behaviour is reflected in the QPC output [Fig. Etc)]: 
when both electrons are localized on dot '1' ('2'), the 
current is smaller (larger) than when the system is in a 
singly occupied state. This is reminiscent of a random 



FIG. 3: Another sample evolution. In this measurement 
run, the state of the system becomes localized in the triplet 
subspace. Parameters as in Fig. 



telegraph switching (RTS) signal. At the first jump, the 
system becomes localized in the singlet subspace [Fig. 
Etb)]. This is because the doubly occupied state |di)(di| 
is a spin-singlet state. Doubly occupied triplet states are 
energetically forbidden, and therefore when a transition 
into the doubly occupied subspace is observed, one can 
infer that the total spin of the system is a spin-singlet. 
In Fig. no jump into the doubly occupied subspace is 
observed, and the detector output current remains con- 
stant (up to shot noise). This leads to a gradual increase 
in the observer's confidence that the CQD system is in 
a spin-triplet state [Fig. OJb)]. In both cases, the final 
state of the system is localized in either the singlet or 
triplet subspace, and a strong measurement of the total 
spin of the system in the singlet-triplet basis is obtained. 

Detector power spectra. The possible detector out- 
puts can be characterized by the power spectrum of 
the current, I(t), through the QPC. The power spec- 
trum is given by S(cu) = 2 f_ dTG(r)e _I "', where 
G{t) = E[I(t+T)I(t)]-E[I(t+T)]E[I(t)], and£LJ de- 
notes the classical expectation. Following Refs. |l3l IbH 
EH, we have G(r) = e 2 [tr{£ n ,„, J[c n ]e CT J[c' n ] Poo } - 

t r {I]„i7MPoo} 2 + tr-E„^MPco}£(T)], where p^ is 
a steady state solution of the UME. Spectra correspond- 
ing to the different measurement outcomes can be found 
by evaluating G(r) for different poo, corresponding to 
steady states localized in the singlet or triplet subspaces. 

For a singlet state outcome, the steady state is ps,oo — 
[(V + U+ J)\s )(s \ + (V -U - J)\ Sl )( Sl \ + (V -u - 
J)\ s 2)(s2\]/(3V — U — J). In the parameter regime of 
interest (cj, J, v 2 V <C J + U < V), the power spectrum is 
well approximated by 



S s {io) =2el + 



mi 2 J 2 v 2 V 



v-u-.i 

3V-U-J 



(w 2 - J 2 ) 2 +4w 2 ( I / 2 V) 2 ' 



(10) 



where I — e 2 (T + v) 2 V is the average current through 
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FIG. 4: Normalized steady state power spectra for the QPC 
detector current for both singlet and triplet steady states. 
Parameters as in Fig. [5] 

the detector, and SI = e 2 (2T v)V . In the limit J < is 2 V, 
Eq. (|10fl corresponds to a RTS process with switching rate 
^rts = J 2 /{^ v2 y) which is the rate at which the system 
jumps between the states \df) and \da). For a triplet state 
outcome, the steady state is an arbitrary state in the 
(singly occupied) triplet subspace. The power spectrum 
only contains a shot noise component, and is given by 
St(<jj) — 2eJ. Power spectra for both singlet and triplet 
outcomes are illustrated in Fig. 

Measurement time and fidelity. The total time re- 
quired to perform a single run of the singlet-triplet mea- 
surement scheme, i mea s has two components, i moa s = 
trel + tdet) where t le \ is the time taken for a singlet 
state to relax into the doubly occupied subspace, and 
tdet = £rts ~ Av 2 V/J 2 is the time required to determine 
the presence or absence of the RTS signal. By analysing 
Eq. we find tf c \ « v 2 J{V - U)/U. i mcas can be 

minimized by varying J(t) over the course of the mea- 
surement, such that different values of J(t) are obtained 
during the 'relaxation' and 'detection' phases of the mea- 
surement. In order to minimize the t IC \, J rc i should be 
as large as possible, but to obtain a strong RTS signal, 
we require Jdet ^ v 2 V . We now evaluate i re i and t<j e t for 
two possible implementations 0, 0] . For P donor qubits 
in Si 0, we take U = 43.8 meV, V = 45.4 meV [23|, 
v 2 = 6 x 10~ 5 (corresponding to a change in conduc- 
tance of 5 % for each electron added to dot T', assuming 
T = 0.3), J re i = 1.0 meV and J do t = 0.8 v 2 V = 2.2 
which gives t re i = 300 ns and tdet = 1.51 ns. 
For quantum dot qubits in GaAs Q, we take U = 1 
meV, V = 2 meV, v 2 =6x 10~ 5 , J rc i = 0.1 meV and 
^dct = 0.8 v 2 V = 0.11 /xeV, giving t Te \ = 110 ns and 
*det = 34.4 ns. 

The fidelity of the measurement is approximately F ~ 
1 ^ ^measAmix) where t mlx is the mixing time for un- 
wanted transitions between the singlet and triplet sub- 
spaces. Owing to the form of the detector-dot interaction 
in Eq.||2J), the detector back-action does not induce such 



mixing transitions, and therefore the dominant contribu- 
tion to t m i x is due to interactions with the environment. 
Although a lower bound of 200 [is for the singlet-triplet 
relaxation time in a single quantum dot has been mea- 
sured , presently no data exists for the two electron-two 
dot case. Therefore we take t mlK ~ minpi, T2], where 
T\ and Ti are the single spin relaxation and dephasing 
times, respectively (both processes can lead to transitions 
between the singlet and triplet subspaces). For P donors 
in Si, T\ can be hours an d ?2 = 50 ms was recently 
measured gi ym g F ~ 1 — 5 x 10 -6 . For GaAs qubits, 
Ti = 0.8 ms has been measured 0], and it is believed that 
T 2 ~Ti [13; thus F ~ 1 - 10" 4 . These results are par- 
ticularly promising for the implementation of single shot 
singlet-triplet measurements with the fidelity required for 
fault tolerant quantum computation |2l| . 
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